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Abstract. We generalize the mode-coupling theory of supercooled fluids to systems under stationary shear flow. Our starting 
point is the generalized fluctuating hydrodynamic equations with a convection term. The method is applied to a two 
dimensional colloidal suspension. The shear rate dependence of the intermediate scattering function and shear viscosity is 
analyzed. The results show a drastic reduction of the structural relaxation time due to shear and strong shear thinning behavior 
of the viscosity which are in qualitative agreement with recent simulations. The microscopic theory with minimal assumptions 
can explain the behavior far beyond the linear response regime. 



Many complex fluids such as suspensions, polymer so- 
lutions, and granular fluids exhibit very diverse rheologi- 
cal behavior. Shear thinning is among the most-known 
phenomena. Recently, it was found by experiments [1] 
and simulations [2] that supercooled liquids near the 
glass-transition also show strong shear thinning behavior 
They have observed that near the transition temperature, 
the structural relaxation time and the shear viscosity both 
decrease as 7^^, where 7 is the shear rate and v is an ex- 
ponent which is less than but close to 1 . For such systems 
driven far from equilibrium, the nonequilibrium param- 
eter 7 is not a small perturbation parameter but plays a 
role more like an intensive parameter which character- 
izes the "thermodynamic state" of the system[3]. Such 
rheological behavior is interesting in its own right, but 
understanding the dynamics of supercooled liquids in a 
nonequilibrium state is more important because it has 
possibility to shed light on an another typical and perhaps 
more important nonequilibrium problem, non-stationary 
aging. Aging is the slow relaxation after a sudden quench 
of temperature below the glass transition temperature. In 
this case, the waiting time plays a similar role to (the 
inverse of) the shear rate. Aging behavior has been ex- 
tensively studied for spin glasses (see Ref.[4] and ref- 
erences therein). The relationship between aging and a 
system driven far away from the equilibrium was consid- 
ered using a schematic model based on the exactly solv- 
able p-spin spin glass by Berthier, Barrat and Kurchan[5] 
and its validity was tested numerically for supercooled 
liquids [6]. There are attempts to study the aging of struc- 
tural glasses theoretically [7] but it has not been analyzed 



and compared with the simulation results [8]. 

In this paper, we investigate the dynamics of super- 
cooled liquids under shear theoretically, by extending 
the standard mode-coupling theory (MCT). We start with 
generalized fluctuating hydrodynamic equations with a 
convection term. Using several approximations, we ob- 
tain a closed nonlinear equation for the intermediate scat- 
tering function for the sheared system. The theory is ap- 
plicable to both normal liquids and colloidal suspensions 
in the absence of hydrodynamic interactions. Numerical 
results will be presented only for the colloidal suspen- 
sions, but generalization to liquids are straightforward. 
Some of the preliminary results have already been pub- 
lished in Ref.[9]. 

We shall consider a two dimensional colloidal suspen- 
sion under a stationary simple shear flow given by 



vo(r) =r-r= (7y,0), 



(1) 



where {r)aii = 75a.x5/3y is the velocity gradient ma- 
trix. The hydrodynamic fluctuations for density p(r,f) 
and the velocity field v(r,f) obey the following set of 
Langevin equations[10]. 



dt 



8^ 

mV ■ (pvv) = -pV- CoP(v - vo) +f«, 

dp 

(2) 



where 1^0 is the collective friction coefficient for colloidal 
particles. fff(r,r) is the random force which satisfies the 



fluctuation-dissipation theorem of the second kind (2nd 
FDT); 

{fR{r,t)fR{r',t'))o = 2k,Tp{r)^o8{r-r')5{t-t') (3) 

for t > t', where (• • • )o is an average over the conditional 
probabihty for a fixed value of p (r) at f = f '. Note that the 
random force depends on the density and thus the noise 
is multiplicative. We assumed that the 2nd FDT holds 
even in nonequilibrium state since the correlation of the 
random forces are short-ranged and short-lived, and thus 
the effect of the shear is expected to be negUgible. The 
friction term is specific for the colloidal case. In the 
case of hquids, it should be replaced by a stress term 
which is proportional to the gradient of the velocity field 
multiplied by the shear viscosity. Both cases, however, 
lead to the same dynamical behavior on long time scales. 
The first term in the right hand side of the equation for the 
momentum is the pressure term and ^ is the total free 
energy in a stationary state. Here we assume that the free 
energy is well approximated by that of the equilibrium 
form and is given by a well-known expression; 



/3^~ |drp(r){lnp(r)/po-l} 

-^y^dri Jdr2 5p(ri)c(ri2)5p(r2), 



(4) 



where /3 = 1 /hgT and c{r) is the direct correlation func- 
tion. Under shear, it is expected that c(r) will be distorted 
and should be replaced by a nonequilibrium, steady state 
form CNE(r), which is an anisotropic function of r. It 
is, however, natural to expect that this distortion is very 
small in the molecular length scale, which plays the most 
important role in the slowing down of the structural re- 
laxation near the glass transition. We confirmed this by 
numerical simulation[l 1]. By linearizing eq.(2) around 
the stationary state as p = po + 5p and v = vo + 5v, 
where po is the average density, we obtain the following 
equations. 



1^ - k-r - ^ ) 5pk(/) = -ikMt), 



|-k.r.|^+k.r.k)/k(0 



ik 



-Spk{t) 



m^S{k) 

--^ / ik-qc(^)5pk-q(05pq(0-— A(0+f«k(0> 
ffip Jq in 

(5) 

where c{q) is the Fourier transform of c(r), k = k/|k|, 
Jk{t) = Pok • 5vk(f ) is the longitudinal momentum fluc- 
tuation, and /q = / dq/(27r)^. Note that our approximate 



equation does not contain coupling to transverse momen- 
tum fluctuations even in the presence of shear. 

In order to construct equations for the appropriate cor- 
relations from the above expressions, an approximate 
symmetry is necessary. In the presence of shear, transla- 
tional invariance is violated. In other words, correlations 
of arbitrary fluctuations, f{r,t) and g{r,t), do not satisfy 
(/(r,f)g(r',0)) ^ (/(r-r',0^(0,0)). Instead, it has die 
following symmetry [12]; 

(/(r,Og(r',0)> = (/(r-r' (0,0^(0,0)), (6) 

or in wavevector space 

(/k(0^k'(0)) = (/k(0^kw(0)) X 5kW,k', (7) 

where we defined the time-dependent position and wave 
vector by r (r ) = exp [Ft] • r = r -h jtyex, where ex is an unit 
vector oriented along the x-axis and k(f ) = exp['n] • k = 
k -I- flkx^y, where T denotes the transpose of F and 
5k,k' = (2;r)^y^'5(k — k') for a system of volume V. 
Using this approximation, it is sttaightforward to con- 
struct the mode-coupUng equations for the appropriate 
correlation functions. We shall derive the equation for 
the intermediate scattering function defined by f (k,/) = 
^~^(5Pk(-;)(^)5Pk(0)), where A'^ is the total number of 
the particles in the system. Note that the wave vector in 
dpk{t) is now replaced by a time-dependent one k(— f). 

Eq.(5) has a quadratic nonlinear term in 5pk(f )■ This 
term can be renormalized to give a generalized fric- 
tion coefficient or the memory kernel following the 
standard procedure of derivation of the mode-coupUng 
equation[13]. To the lowest order in the loop expansions, 
we obtains the equation for the velocity-density correla- 
tionC(k,0=iV-i(4(_,)(;K(0)); 



dC(k,0 



■k(-?)-r-k(-OC(k,0 : 



ik{-t)F{k,t) 



m^S{k{-t)) 

-- fdt' fdk.' ak(-t),k',t-t')C(k',t'). 
m Jo J 

(8) 

Note that in the above equation, the differential op- 
erator k-F-d/dk disappears and k is replaced by 
k(/). C(k,k',/) is the generalized friction coefficient. 
C(k,k',?) is given by the sum of the bare friction coeffi- 
cient and the mode-couphng term as 

C{k,k',t) = Co X 25{t) + 5C(k,05k(0,k' (9) 
with the mode-coupling contribution given by 



5C(k,O=0_^r(k,q)r(k(f),q(O) 
xF{k{t)-q{t),t)F{q{t),t), 



(10) 



where l^(k, q) is the vertex function given by 

r (k, q) = k • {qc(q) + (k - q) c(k - q)} . (11) 

Note that, in the derivation of eq.(lO), we have assumed 
the fluctuation-dissipation theorem of the first kind (1st 
FDT), which relates the response function to the correla- 
tion function. In order to derive the renormalized friction 
coefficient, propagators (the response function to the ran- 
dom forces) as well as correlation functions [13] are natu- 
rally introduced. The 1st FDT makes it possible to elim- 
inate the propagators in favour of the correlation func- 
tions. In the overdamped limit, we may neglect the time- 
derivative of C(k,/). The second term on the left hand 
side is also neglected if Peclet number Pe= ya'^/D^ {a 
is the diameter of the particle and Do = k^T/Q) is the 
diffusion coefficient) is small. Therefore, combining the 
first equation of eq.(5), one may eliminate C(k,f) from 
the above equations and arrive at the closed equation for 
F(k,f); 



dF{k,t) Dok{-tf 



S{k{-t)) 



(12) 



where 



MiKt)='-^^^fnK^)nm,m ^^^^ 

xF{k{t)-qit).t)F{q{t)j). 

In the absence of the shear, they reduce to the conven- 
tional mode-couphng equation[14]. 

Solving eqs.(12) and (13) numerically is more de- 
manding than the equations in the equihbrium state be- 
cause the wavevectors are distorted by shear and the sys- 
tem is not isotropic. We have considered a two dimen- 
sional colloidal suspension consisting of hard disks. For 
the static correlation function c{k) and S{k), the analytic 
expressions derived by Baus et aJ.[15] were used. We 
have divided the two dimensional wavevector space into 
Nk grids for each direction. The cut-off wavevector was 
chosen to be kcO = IOtt. For the self-consistent calcu- 
lation of the mode-coupling equation, we have used the 
algorithm developed by Fuchs et aJ. [ 1 6] . In the following 
results, we have used the grid number Nk = 55, for which 
the ergodic-nonergodic transition occurs at the volume 
fraction (pc = Tta^po/A = 0.7665. Apparently, Nj^ is not 
large enough to give the right transition density which is 
10% smaller. However, the quahtative behaviors does not 
change by increasing the grid number. 

In Figure 1, we show the behavior of F(k,f) for = 
X (1 - lO""^) for various shear rate Pe=10"'° to 10. 
The wavevectors were chosen to be kc7 = (0,3) and 
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FIGURE 1. F{\i,t)/S{k) for ktr = (0,3), (3,0) for various 
shear rates y. From the right to the left, Pe= 10"'*', 10~', 10~^, 
10"^, 10"', and 10. The time t is scaled by o^/Dq. 



(3,0), parallel and perpendicular to the shear flow, re- 
spectively. For the smaller shear rate, Pe < 10"'", we did 
not see any shear effect. For higher shear rates, we have 
observed the drastic reduction of the relaxation time due 
to shear. This is similar to the behavior reported in re- 
cent molecular dynamics simulations for a binary soft- 
core liquid[2]. 

The two lines for a fixed shear magnitude but for 
wavevectors along different directions collapsed onto 
each other and we do not see a noticeable difference. It 
is surprising that, though the perturbation of shear flow 
is highly anisotropic, the dynamics of fluctuations are al- 
most isotropic. This fact is also observed in the simula- 
tion for a binary liquid[l 1]. The reason for the isotropic 
nature is understood as follows. The shear flow perturbs 
and randomizes the phase of coupUng between different 
modes. This perturbation dissipates the cage that tran- 
siently immobilizes particles. Mathematically, this is re- 
flected through the time dependence of the vertex. This 
"phase randomization" occurs irrespective of the direc- 
tion of the wavevector, which results in the isotropic be- 
havior of relaxation. This mechanism is very different 
from that of many complex fluids and dynamic critical 
phenomena under shear, in which the faster relaxation 
occurs mainly due to the distortion of the structures at 
small wavevectors which are stretched out by the shear 
flow and pushed to larger wavevectors where faster re- 
laxation occurs. 

The shear viscosity, t], is easily evaluated by modify- 
ing the Green- Kubo formula for the sheared system[17]; 



'?(7) = '7o 



kMt) dSjk) dS{k{t)) 

kS2{k)S^{k{t)) dky dky{t) 



F\k{t),t), 
(14) 



where tjo is the viscosity of the solvent alone. The inte- 
gral can be implemented for the set of F{k,t) evaluated 
using eq.(12). Shear rate dependence of the reduced vis- 
cosity r\R{f) = {77 (7) - 77o}/tJo is plotted in Figure 2 for 
various densities around ^c- The strong non-Newtonian 
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FIGURE 2. The reduced viscosity is plotted for Pe = 
ya^ /Dq for various densities. From above to bottom; = 
0.766549, 0.7665, 0.766453, 0.7664, 0.766, and 0.756. The 
highest density is 4 x 10~^ % larger than 0^. 

behavior is observed at high shear rate and large densi- 
ties, which is again in qualitative agreement with the sim- 
ulation results for liquids. Slightly above the plastic 
behavior which implies the presence of the yield stress is 
also observed. The shear thinning exponent is extracted 
from this plot between 10^ < Pe < 1 and we obtained 
f]R{y) °^ with V ~ 0.99. For the larger shear rate, Pe 
> 1, the exponent becomes smaller. In this regime, it is 
expected that other mechanism such as the distortion of 
structure c{k) and S{k) by shear becomes important. 

The mode-coupling theory developed in this paper is 
far from complete. The most crucial approximation is 
the use of the 1st FDT, which was employed when we 
close the equation in terms of the correlation functions 
alone. It is already known that the 1st FDT is violated 
for supercooled systems under shear as well as during 
aging [6]. Without the 1st FDT, one has to solve simulta- 
neously the set of mode-coupling equations for the prop- 
agator and correlation function, which couple each other 
through the memory kernels. Research in this direction is 
under way. Another important approximation was to ne- 
glect the small distortion of the structure, c{k) and S{k), 
due to shear. The construction of the equation for such 
an equal-time correlation functions might be more sub- 
tle and should be considered in future. It is surprising, 
however, that despite of these approximations, the theory 
reproduces the major features which was seen in simu- 
lations; the drastic reduction of the relaxation time, the 
isotropic nature of the dynamics, and plastic-hke strong 
shear thiiming. 



An analogous effort has been made by Fuchs and 
Gates [18]. They have derived a mode-coupling expres- 
sion for F{k,t) using a projection operator for the Smolu- 
chowski equation for A^-particle colloidal suspensions. 
They have observed similar shear thinning behavior for 
approximated expressions, for an isotropic model, where 
the anisotropy hidden in the equations are neglected. 

The details of analytical and numerical calcula- 
tions for results given in the present paper are given 
elsewhere[ll]. 
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